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ABSTRACT 

It has been suggested that Narrow-Line Seyfert 1 (NLSl) galaxies are evolutionarily young objects, 
powered by the accretion of gas onto central black holes that are significantly lower in mass than those 
found in typical broad- line Seyferts. We explore this hypothesis through the analysis of high-spatial 
resolution, ncar-IR imaging data obtained in J and K' for a sample of 11 NLSls. Surface bright- 
ness profiles are separated into their constituent components using two-dimensional decomposition 
techniques. By employing the correlation between black- hole mass and host galaxy bulge luminosity, 
calibrated for near-IR wavelengths using 2MASS data, we determine the mean black-hole mass for 
our sample to be, in solar units, (log(Ms_ff)) = 7.9. Using the correlation between the size of the 
broad-line region and the monochromatic continuum luminosity, we obtain black-hole mass estimates 
under the assumption that the emission-line gas is in virial equilibrium. The mean black-hole mass 
derived from this relation is (log(MB//)) = 6.4. It is found that the estimates obtained from the 
black-hole mass-bulge luminosity relation are systematically one full order of magnitude larger than 
those derived from the black-hole mass-broad-line region radius relation. We explore possible causes 
for this discrepancy in Mbh estimates and the ramifications for our understanding of the role played 
by NLSls in AGN evolution. 

Because numerical simulations constrain the start of the AGN duty cycle to a time shortly after a 
significant gravitational interaction, we examine the morphology and near-IR bulge colors of the NLSl 
sample for evidence of recent encounters. The mean bulge color is found to be {{J — Kg)) = -1-1.85 
± 0.58, which is redder than that of both a matched sample of non-active galaxies and published 
estimates for broad-line Seyferts. The source of the unusual bulge colors may be an excess of flux, 
peaking at around 2.2 /im, that has been detected near the centers of some NLSls such as Mrk 1239. 
No evidence is found for light asymmetries or an extra stellar component that would indicate NLSls 
are young objects. Finally, we postulate that there may be some interesting lines of circumstantial 
evidence suggesting that secular processes may be relevant in the evolution of NLSls. 
Subject headings: galaxies: Seyfert — galaxies: evolution — galaxies: photometry — infrared: galaxies 
— techniques: photometric 



1. INTRODUCTION 

The connection between the evolution of galaxies and 
that of the supermassive black holes that are known to 
lurk at their centers is a key issue in contemporary extra- 
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galactic studies. It is clear that Active Galactic Nuclei 
(AGN) represent an important, transient phase in t he de- 
velop ment of galaxies at all redshifts (Hachnclt fc R.eetj 
ri99l . The source of intense, non-thermal emission gen- 
erated in the cores of AGN is widely accepted to be the 
accretion of material onto otherwise d ormant black holes 
(|R,eeslll984l IPeterson fc Wandelll200nD . However we may 
not yet fully appreciate the diversity of mechanisms that 
create the non-equilibrium conditions leading to accre- 
tion. Significant changes in gravitational potential can be 
produced by interactions between galaxies of comparable 
mass or by the consumption of dwarf satellite galaxies 
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l)Osterbrocklll993j) . It also may be possible to trigger ac- 
cretion through secular processes within the host galaxy 
such as bar-driven tidal instabilities, generating episodic 
bursts of nuclear activity (Kormcndv & Kcnnicutt 2004). 

An evolutionary link between black holes and their 
host galaxies is strongly suggested on the basis of cor- 
relations between the black-hole mass, Mbh, and var- 
ious galaxy parameters. A good correlation has been 
found between Mbh and the central stellar velocity dis- 
persio n, g.,, of the sp heroid compo nent of the host galaxy 
(e.g.. iFerrarese fc Merritt 200Qe iGebhardt eraP 1200(1 
iGreene fc Holl2005ajl . The intrinsic scatter in this re- 
lation ma y be more significan t at low black-hole masses, 
however ijBender et al.ll2005l) . It has also been shown 
that Mbh is a constant fraction of the mass of the host 
spheroid (or, in the case of ellipticals, the total mass of 
the galaxy), with Mbh I M^r^h — 0.0015 (Magorria n et al.l 
ITqM IMarconi fc HunUMB IHaring fc Rix..2004l) . Un- 
der the assumption of a constant mass-to-light (M/L) 
ratio, this relation is frequently expressed in terms of the 
host bulge luminosity (e.g. Lao r 2001: Mc Lurc & Dunlonl 
120011) . While early studies found the scatter in this mass- 
luminosity relation to be substantial, recent progress in 
the amount and quality of available data has led to 
an irnprovement in the precision of the best-fit parame- 
ters llMarconi fc Hunt 2003; Dong fc Dc Robertis 200^. 
IMarconi fc Hund lf200.'^l also demonstrated the wave- 
length dependence of the scatter in the relation, with 
near-IR observations tracing the mass distribution bet- 
ter than optical data, generating a tighter correlation. 

There is a general consensus in the literature that AGN 
exhibit these correlations independent of distance or lu- 
minosity, suggesting that they are a natural product of 
galax y-black hole co-evohiti on llKormendv fc Richstoii^ 
119951 iBarth. Greene, fc Hoi 120051^ Yet there remains 
much debate as to how Narrow-Line Seyfert 1 galax- 
ies, which constitute approx imately 15% of the total low- 
redshift Seyfert population l|Williams. Pogge. fc MathuJ 
|2002), fit into this picture. NLSls are a subset of AGN 
exhibiting an overall spectral signature consistent with 
a Type 1 Seyfert, but with remarkably narrow permit- 
ted emission lines. More specifically, NLSls are classi- 
fied based on an observed broad component of H/? emis- 
sion with FWHM < 2000 km s^^, strong Feii emission, 
relatively weak [O iii] , and steeper soft x-ray continuum 
slopes than broad- line Seyfert Is (jOstcrbrock & Pogg9 
119851 iBoller. Brandt, fc Fink 1996). _Througha Princi- 
pal Component analvsis. Boroson fc GreenI l)1992t) found 
that NLSls dominate the extreme end of the "Eigen- 
vector 1" anti-correlation between the Feii and [Oiii] 
emission-line strengths. This relation is believed to be 
driven by the ratio of the luminosity of t he source rela- 
tive to its Eddington luminosity, L/LEdd ijSulentic et al.l 
I2OOO). Clearly, NLSls occupy an important portion of 
AGN parameter space. 

To account for the extreme properties of NLSl galax- 
ies, it has been conjectured that their central engines 
are black holes of relatively modest ma s s (ty pically of 
order 10*^ M©) l)Boller. Brandt, fc FinkI [1996 : Mathur 
l2000(l . This hypothesis can explain the narrow emis- 
sion lines observed in these objects, under the assump- 
tion that the broad H/3-emitting region is gravitation- 
ally dominated by the black hole. If the black holes 
in NLSls are statistically smaller than in broad-line 



Seyfert galaxies, elevated accretion rates would be re- 
quired to a ccount for the relatively normal o bserved lu- 
minosities l|Pounds. Done, fc Osborii^ Il995j) . Further- 
more, it has been suggested that these objects repre- 
sent an early stage in AGN evolution, during which the 
black-hole mass increase s rapidly in a period of intense 
accretion ()Mathuill2000|) . As the black hole grows m 
this scenario, the accretion rate gradually decreases and 
the NLSl matures into a "normal" broad-line Seyfert 
ijCrupe fc M athur '2 0041) . 

It remains unclear whether the aforementioned cor- 
relations with Mbh can be applied to NLSls. Sev- 
eral authors suggest that NLSls fall below the mass- 
luminosity correlation of other galaxies, resulting in 
smaller estimates of Mbh for a given bulge luminosity 
(Mathur. Kuraszkicwic z. fc Czernv 2001; Wandcl 2002). 
It has also been suggested that intense star formation ac- 
tivity could lead to overestimates of the bulge luminos- 
ity or, equivalently , underestimates of the {M/L) ratio 
(|Nelson ct al. "2Q0^. The vahdity of the Mbh-cTs corre- 
lation for NLSls has also been questioned. It has been 
demonstrated bv lBarth. Greene, fc Hoi l|2005D that AGN 
with small black holes follow the same relat ion as galax- 
ies wit h m ore typical values fo r Mbh- Both I Wang fc Ld 
lj200H) and lBotte et alJ lj2004 t) estimate the velocity dis- 
persion for NLSl galaxy samples from measurements 
of [Oiii], and find that the Mbh-Cs correlation is the 
same as for broad- line Seyferts. Conversely, a depar- 
ture from th e standard Mbh-cTs relation for NL S ls wa s 
reported by 'Mathur. Kuraszkicwicz. fc Czernv' ('2001), 
GruDC & Mathur (2004), and Bian & Zhao (2004). Each 
of these studies used [Om] as a surrogate for the ve- 
locity dispersion. It is interesting to note that different 
conclusions are reached in these studies even though the 
samples are similar. 

A major thrust of this study is to compare black-hole 
masses derived from the well-known correlation between 
Mbh and the host galaxy luminosity for the NLSl sub- 
set, hereafter the mass- luminosity correlation, with other 
mass estimates. To this end, we examine a sample of 
11 NLSls using high-spatial resolution, near-IR imaging 
data. The luminosities of the bulge components in J and 
K' are measured independently for each target. Black- 
hole masses are obtained by utilizing the mass-luminosity 
relation, calibrated to near-IR wavelengths. This work 
builds on previous studies by applying two-dimensional 
profile decomposition techniques to near-IR data of ex- 
cellent quality, providing an accurate determination of 
the photometric properties of the host bulges. The re- 
sults are compared to values of Mbh derived via the 
correlation between the radius of the broad-line region 
and the optical continuum luminosity, which has been 
calibrated us ing reverberation mapping mass estimates 
l|KasDi et al.l l2000. 200^). This comparison is used to 
examine the role of NLSls in the evolution of AGN. 

In § 2 we outline the sample selection and describe 
the observations. The data reduction process, including 
the decomposition procedure, is included in § 3. The re- 
sults and a discussion of their implications for the NLSl 
paradigm are presented in § 4 and § 5, respectively. Con- 
clusions are provided in § 6. Throughout this work, we 
adopt cosmological parameters of JIa — 0.7, = 0.3 
and a Hubble constant of Ho = 75 km s~^ Mpc~^. 
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2. OBSERVATIONS 

To examine the possible role of NLSl galaxies in the 
evo lution of nuclear activity, a sa mple was selected from 
the|V eron- Cettv fc Vero^ l|200H) catalog, Quasars and 
Active Galactic Nuclei, 10th Edition. Targets were 
chosen to be low-redshift (z < 0.05) and moderately 
bright (total K' < +12.0) in order to achieve the excep- 
tional spatial resolution and signal-to-noise ratios nec- 
essary for accur ate profile decomposit ions. Selecting 
targets from the iVeron-Cettv fc VeronI catalog also en- 
sured that a reasonable amount of complementary ob- 
servations wou ld be available, including optical spec - 
troscopic data ijVeron- Cettv. Veron. fc rxoncalve!^l2?ffll . 
A total of 150 object s are classified as NLSls in the 
iVeron-Cettv fc VeronI catalog. Of these, approximately 
37 matching our selection criteria are visible in the 
Northern sky. We imaged 11 of these NLSls in De- 
cember 2002 through the broad-band J and K' filters 
at the 3.6-meter CFHT using KIR, the 1024 x 1024 
HgCdTe array and near-IR camera, mounted at the f/20 
output focus of the PUEO Adaptive Optics Bonnette 
l)Rig-aut et al.f lQQS'). The focal plane scale of this system 
is 0.0348"/pixel, providing a total field of view of 36" x 
36". For all of the target observations, the galaxy nu- 
cleus was sufficiently bright to be used as the wavefront 
reference source. The photometric standard stars were 
also observed in the same way, acting as their own wave- 
front reference sources. This greatly facilitated the reg- 
istration and co-addition of all dithered sequences, ren- 
dering other corrections unnecessary. Strehl ratios for 
the galaxy observations were estimated to be between 
0.2 and 0.4, while for the standard st ars the Strehl ratio 
was found to be approximately 0.5 l|Beuzit fc Hainaxil 
[1998; Rigaut et al. 1998). A PSF FWHM of 0.2" was 
consistently measured from the standard star exposures, 
providing significant spatial over-sampling of the data. 
The key photometric properties of the galaxies, includ- 
ing morphological type, optical semi-major axis, and K 
— +20 mag arcsec"^ isophotal radius, are provided in 
Table nj We note that the if-band isophotal radii of the 
targets are between 6.4" and 26.1". 

An exposure sequence for each object consisted of six 
spatially-dithered images, including two off-target expo- 
sures that were subsequently used to carefully monitor 
the near-IR sky and to provide accurate background sub- 
traction. The data were acquired using the sequence T, 
T, B, T, T, B, with T representing an on-target ex- 
posure and B an off-target background exposure. The 
off-target exposures were dithered by between 1.6' and 
2.3', depending on the size of the galaxy, to the East 
and West sides of the on-target field of view in order to 
avoid contamination by the extended disks of the galax- 
ies. This procedure ensured that the on-target images 
were bracketed, both spatially and temporally, by the 
background frames. The on-target exposures were typi- 
cally dithered by approximately 16" in a four-point box 
pattern that placed the galaxy nucleus in each of the de- 
tector's quadrants. Total on-target exposure times var- 
ied between 80 and 1920 seconds, chosen to maximize 
the signal while ensuring that the detector response re- 
mained linear. Equal exposure times were used for all 
images in each sequence. The mean airmass for each 
data set was typically between 1.0 and 1.2. A summary 



of the observation log is provided in Table |21 
3. DATA REDUCTION 

The data were reduced using standard IRAF proce- 
dures, including the construction and application of bad 
pixel masks, the removal of dark current, and a flat-field 
correction (including any large-scale gradients). A cru- 
cial step in the reduction and analysis of near-IR images 
is the accurate subtraction of the sky background. The 
characterization of the background level and pattern is 
particularly challenging, as it has been shown that the 
mean level at the CFHT site can vary by 0.7% and 0.5% 
in J and K\ respectively , on time scales of one minute 
ijVaduvescu fc McCaU 12004) . Failure to take this into 
consideration could potentially lead to photometric er- 
rors when imaging extended sources with lengthy inte- 
gration times. Because of the relatively large angular 
size of the targets compared with the detector's field of 
view, it was impractical to determine the background 
level using adjacent on-target exposures. Instead, accu- 
rate background subtraction was achieved by carefully 
monitoring the two off-target exposures from each data 
set. An inspection of the median and mode for each 
pair of off-target frames indicated that the background 
level and underlying patterns were constant to within 
1%, allowing us to use the combined off-target frames 
for background subtraction. 

We then examined the on-target exposures, both be- 
fore and after background subtraction, to search for any 
systematic uncertainties introduced by this procedure. 
As the on-target dither pattern situated the galaxy nu- 
cleus sequentially in each of the detector's four quad- 
rants, we monitored the image statistics of the opposite 
quadrant carefully. Variations in the median count rate 
among the image corners were always within underlying 
statistical uncertainties. The residual ffux levels of the 
corners of the background-subtracted images were con- 
sistently near zero, with no residual gradients present. 
This gave us confidence that the background had been 
successfully and appropriately removed in all target ex- 
posures. 

The dithered on-target exposures were combined into 
56" X 56" mosaic images using the IRAF package ir- 
red.irmosaic. For three galaxies, multiple image sets were 
obtained in K', allowing independent mosaics to be gen- 
erated. This served as a useful check on the consistency 
of the data reduction and analysis processes. 

3.1. Bulge-Disk Decomposition 

The surface brightness profiles of the reduced mo- 
saic images were then examined to determine the rela- 
tive flux associated with the bulge component. To ac- 
complish this, the two-dimensional proflle fitting algo- 
rithm GALFIT v2.0 was employed. GALFIT allows the 
user to model the light profile of a target using a suite 
of analytic functions, including Gaussian, exponential, 
Sersic/de Vaucouleurs, and Nuker profiles. The optimal 
fit is determined, based on the user's selection of input 
parameters, by minimizing the residual between orig- 
inal and model profiles. Details of the minimization rou- 
tine are provided in Pcne et al. (2002). 

While GALFIT is capable of fitting a PSF-convolvcd 
Gaussian to account for a point source, in the interest of 
limiting the number of free parameters in the model and 
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to best ensure the parameters are physically realistic, an 
alternative method was used to subtract the contribution 
of the active nucleus to the integrated flux. This shift- 
scale-subtract (sss) routine involves using a high signal- 
to-noise ratio, isolated point source such as a photo- 
metric standard sta r observed under similar conditions 
to m odel the AGN l|Virani. De Robertis. fc VanDalfsenI 
|200C1|) . For each target, a standard star observed close in 
time was selected as a template point source. The stars 
were imaged with the AOB activated, using the same 
dither pattern as employed for the target observations. 
The spatial and temporal stability of the standard stars 
was monitored by examining the reduced exposures. It 
was found that the stellar PSFs were consistent to within 
0.02 magnitudes, and the FWHM varied by less than 
0.05". To subtract the nuclear component, the star was 
aligned to the center of the galaxy as estimated from el- 
liptical isophote fitting via the ellipse package. Accura- 
cies of ±0.1 pixel are achievable for the spatial alignment 
by monitoring asymmetries in the light profile of the dif- 
ference image. The scaling of the nuclear PSF was judged 
to be that which produced a nearly flat profile in the cen- 
tral 0.4". The average contribution of the nuclei to the 
total galaxy flux was determined to be approximately 
15% in J and 9% in K' . It was thus possible to assess 
the nuclear contribution to better than 1% of the to- 
tal galaxy luminosity. Appropriate nuclear scale factors 
were determined independently for the J- and K'-hand 
images. 

The nucleus-subtracted galaxy images were then mod- 
eled as the sum of a bulge and, when necessary, disk 
component using GALF IT. Galaxy dis ks were fit using 
an exponential profile (jFreemanI 1 19701) . The analytic 
form of galaxy bulges is less clear, and has been dis- 
cussed extensively in the literature. Most early stud- 
ies represented bulge light as a de Vaucoulcurs (1948) 
7-1/4 pj-ofile, similar to that used for elliptical galax- 
ies. However, for many galaxies the profile fails to 
accur ately reproduce the shape of the bulge at all 
radii (jCaon. CapaccioU. k D 'Onofriol Il993j) . This de- 
parture from the r^^^ profile is of greatest concern for 
early-type galaxies, where a significant fraction of the 
total galaxy light comes from the bulge component. 
It has been demonstrated that bulg es ar e better de- 
scribed using the more gener al Sersi 4jl96^ ri/" profile 
ijAndredakis. Peletier. fc BaTcellsii il99,'iD . which reduces 
to a de Vaucouleurs profile for n = 4. Sersic profiles 
were applied to model the bulge light for the NLSl sam- 
ple with the index, n, a free parameter to be determined 
by GALFIT. As all of the targets are early- type galaxies, 
no additional components were required to account for 
features such as spiral arms. 

Initial estimates for the disk scale radii were taken 
from the 2MASS Extended Source Catalog (XSC), re- 
calibrated to the KIR image scale. The XSC images 
have a relatively coarse l"/pixel focal plane scale and 
an integration time of 7.8 seconds. As such, the pub- 
lished scale parameters for these moderately faint ob- 
jects, which have a mean 2MASS Kg = -1-10.8, were 
treated as initial estimates to be determined more ac- 
curately by GALFIT. 

Starting values for the bulge scale radii were then es- 
timated from an examination of disk-subtracted, one- 



dimensional surface brightness profiles. The position an- 
gle, ellipticity and a4 "diskiness parameter" were also 
estimated from elliptical isophote fitting to the nucleus- 
subtracted images. For the Sersic index, an initial esti- 
mate of n = 4 was used, as the targets are generally early- 
type galaxies and it has been demonstrated that there is a 
trend of decreasing index as morph ological type proceeds 
from early to late l)Grahamll2001|) . All of these param- 
eters were allowed to vary within the minimization 
process, however. The luminosity bulge-to-disk ratio was 
estimated from the empirica l relation with Hubbl e stage 
originally presented in Simi en fc dg_Ygi ^mileiir^ l|198fif) 
and derived from if-band data in iGrahamI 1)200 ID . 

For each target, we used GALFIT to first obtain one- 
component fits to the disks of the galaxies. Pixel masks 
aligned to the isophotal centers of the galaxies were used 
to block the bulge components in these initial model fits. 
The results were then used as estimates of the input disk 
parameters in two-component fits, again allowing all pa- 
rameters to vary in order to find the global minimum. 
For these fits, we did not apply pixel masks to any sec- 
tions of the mosaic images and fit the galaxy profiles to 
r = 0. 

It is also necessary that the relevant parameters be 
physically realistic. For three targets, Mrk 359, Mrk 1044 
and MCG08.15.009, it was necessary to constrain one of 
the free parameters in order to avoid convergence to a 
non-physical analytic fit. For these galaxies, the bulge 
scale radius was held fixed to the value provided by the 
2MASS XSC. Residual images were also visually exam- 
ined for any significant features or asymmetries. 

Proper background subtraction is essential to measure 
accurate integrated fluxes for each component of the 
galaxy. To estimate the uncertainties in the fitted pa- 
rameters to reasonable uncertainties in the background 
level, we performed fits to the data with sky levels set 
systematically high and low by one standard deviation 
in the background. The same input parameters used to 
derive the optimal fits were employed in these tests. This 
technique is similar to that applied to one-dimensional 
decompositions by Graham (2001). In some cases, neg- 
ative fluxes introduced by background over-subtraction 
prevented GALFIT from converging on a meaningful so- 
lution. Based on the results of these extensive tests, how- 
ever, it was established that conservative uncertainties in 
the bulge fits resulting from imperfect background sub- 
traction are less than 0.1 magnitudes. 

Separate model fits to the J and K' images were ob- 
tained for each galaxy. Plots of the J-band surface bright- 
ness decompositions are shown in Fierures ITI throuerh II II 
For each target, the relative contribution of the compo- 
nents are presented in log space; residuals to the fit are 
plotted separately in linear space. For the if'-band de- 
compositions, we find that there is a consistent pattern in 
many of the residual profiles on angular scales less than 
1" from the galaxy center, such that a broad negative 
peak is followed by a small upturn. This is illustrated 
in Figure El which presents the iiT'-band profile for Mrk 
1126. It is conceivable that the pattern, which has an 
insignificant effect on total bulge magnitudes, may be 
due to limitations in our nucleus-subtraction routine. In 
order to search for any systematic effects introduced by 
the sss-technique, we performed profile decompositions 
on images in which the point sources were purposely over- 
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subtracted and under-subtracted by 10%. This level of 
uncertainty is greater than could reasonably be expected 
from a careful subtraction of the nucleus. We found in 
these cases that the bulge magnitudes varied by of order 
±0.03 magnitudes. The reason this has no significance on 
total bulge magnitudes is essentially because the angular 
extent of the model bulge profiles are between 2" and 
8". Since nuclei affect only the central 0.2", their in- 
fluence is limited to a fraction of the inner 1% of the 
two-dimensional bulge profiles. 

In employing the sss-technique, it is assumed that the 
PSFs associated with the active nuclei can be adequately 
modeled as scaled versions of the photometric standard 
stars. However, small discrepancies could arise from ef- 
fects such as changes in the airmass of the observations 
or contamination from the AGN host galaxies. Fur- 
thermore, differences in the wavefront reference source 
magnitudes between the AGN and standard star obser- 
vations resulted in small variations in the Strehl ratios. 
This could produce slight differences in the AO-corrected 
PSFs. To further verify that the subtraction of the flux 
associated with the active nucleus did not introduce any 
systematic effects on the bulge fits, we performed three- 
component fits to the data using GALFIT. In addition 
to the bulge and disk components, a nuclear point source 
was included, represented by a Gaussian function. The 
bulge magnitudes and Sersic indices derived from these 
fits were consistent with the results obtained using the 
sss-technique. The quality of the fits, as parameter- 
ized by the reduced values, were also comparable. 
Additionally, two-component (bulge-|-disk) fits were per- 
formed on the non-PSF-subtracted images. The output 
bulge magnitudes were again similar to the results pre- 
sented in Table |3| albeit with slightly higher Sersic in- 
dices and increased residuals. These results are not sur- 
prising, given that the fiux associated with the nucleus 
is of order a few percent of the total galaxy fiux, and is 
confined to the central 1% of the model fits. From these 
tests it was concluded that the sss-technique is robust, 
and does not compromise the model bulge parameters. 

The small but systematic if '-band residuals may also 
reflect a departure of the analytical model from the physi- 
cal structure near the center of the bulge component that 
is resolved by our data. An additional source of emis- 
sion may influence the surface brightness on scales less 
than 1" from the nucleus ([Rodriguez- Ardila & Mazzalay, 
|2006|) . This option is discussed in greater detail in § 5. 
The net effect of these small, systematic residuals on the 
total integrated luminosities of the galaxy components is 
negligible. 

For the three galaxies that were observed twice in 
K' , independent flts provided a useful check of the con- 
sistency of the fitting procedure. Wc found that the 
bulge magnitudes obtained from these images were con- 
sistent within the formal uncertainties. Absolute bulge 
magnitudes were then calculated using redshift infor- 
mation taken from the literature and Galactic extinc- 
tion corrections based on the dust maps presented in 
ISchlegel. Finkbeiner. fc Davi3 1)19981) . The results are 
summarized in Table 

4. BLACK-HOLE MASSES 

Black-hole mass estimates for AGN are frequently ob- 
tained by employing the relation between the size of the 



broad-line region and the monochromatic continuum lu- 
minosity at 5100 A, the mass-radius relation. We will 
compare these results with masses derived from the well 
known correlation between Mbh and the host bulge lu- 
minosity, the mass-luminosity relation, using the results 
obtained from the two-dimensional profile decomposi- 
tions. In the following subsections we discuss these tech- 
niques in greater detail, and present the resulting black- 
hole mass estimates. 

4.1. The Mass-Radius Relation 

Reverberation mapping is widely considered to be a 
reliable method to obtain accurate Mbh estimates in 
Seyfert 1 galaxies, for which direct probing of the black- 
hole sphere of infiuence is not possible due to the intense 
nuclear emission. Correlations with host galaxy parame- 
ters derived from reverberation mapping results are gen- 
erally consistent with those of non-active gal axies with 
secure values of Mbh l|Ferrarese fc Fordll2005|) . However 
direct measurements of the time delay between fiuctua- 
tions in continuum and line emission can be challenging 
and require long-term monitoring of the source. As a 
consequence, estimates of Mbh are more commonly gen- 
erated using secondary, indirect methods such as exploit- 
ing the correlation between the radius of the broad-line 
region and the monoc hromatic continu um luminosity at 
optical wavelengths fKaspi et al.ll200,^ : 

n ^^r /ri « ^ H 0.69±0.05 

Rblr , s rALA(5100A) 

— — = 2.23 ±0.21 ^ ' 
10 It - days ^ ' 
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Based on a sa mple of six obj e cts w ith FWHM(H/3) 
< 2000 km s-\ iPeterson et"aLl l)2000|) suggested that 
this relation can be extended to the NLSl subset. This 
conclusion is supported by the observation that the con- 
tinuum luminosity at 51 00 A is cornparab le in NLSls 
and broad-line Seyferts l|Grupe et alJl200^ . However, 
the low optical variability that is characteristic of NLSls 
makes secure reverberation mapping masses of these ob- 
jects particularly difficult to obtain, as demonstrated by 
the small number that hav e been successfully monitored 
(e.g. iShemmer et"aI]l200lD . 

Optical spectra were published by 

iVeron-Cettv. Veron. fc G oncalvcs (2001) for a large 
number of NLSls, including seven of the objects in our 
sample. The data were obtained with the CARELEC 
spectrograph mounted to the 1.93-meter telescope at the 
Observatoire de Haute-Provence, providing a resolution 
of 200 km s~^ at H/3. The contribution of Feii emission 
was removed using a template spectrum generated from 
observations of the N LSl I Zw 1, according to the 
procedure discussed in iBoroson fc GreenI l)1992D . The 
published spectra have a wavelength range of 4750 to 
5100 A. The paper also includes values for the FWHM 
of the broad component of H/? for each object, taken 
from the literature. We measured the monochromatic 
continuum luminosity at 5100 A from these spectra, 
correcting for Galactic extinction. The results were 
used to estimate the radius of the broad-line region 
according to Assuming the emission-line gas is in 
virial equilibrium with the central black hole, estimates 
of Mbh for each object may then be obtained using the 
following expression: 

MBH^RBLRG-^ikvHSif. (2) 
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The constant A; is a function of the geometry of the broad- 
hne region. Assuming t he hne-emitti i ig gas is dis tributed 
purely isotropically, as iKaspi et alJ l)200(]|) and iWandell 
((2002), k= V3/2. 

For e ach of our targets with spectroscopic da. ta avail- 
able in iVeron-Cettv. Veron. fc Goncalve'3 1)200 Ij) . values 
for the continuum luminosity at 5100 A, the FWHM of 
the broad component of H/?, and the calculated black- 
hole masses are included in Table |5l The mean value of 
Mbh for the seven NLSls is {\og{MBH)) = 6.3, in agree- 
ment with the idea that black holes in NLSls are statisti- 
cally less massive than broad-line Seyferts. The assump- 
tion of a spherically symmetric gas distribution, however, 
is likely to be an oversimplification. A physically more 
meaningful model may i nclude a gas disk in addit ion to 
an isotropic component lIMcLure fc Dunlodl2(ffll . The 
introduction of a disk component in the broad-line region 
would increase the virial mass estimates presented here 
b y a factor of th ree . 

iGreene fc Ho ( 2005b) propose a correlation between 
the radius of the broad-line region and Balmer emission- 
line luminosities. Using a sample of 229 objects take n 
from the Sloan Digital Sky Survey (jYork et al.ll2000j) . 
they obtain an empirical relation between Mbh and the 
H/3 line luminosity of the form: 



Mbh = (3.6±0.2)xl0^ 



Lh(3 
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0.56±0.02 



FWHM 



H/3 



lO^kms"^ 

In formulating this relation, the aut hors derive the mass - 
luminosity relation independent of IKaspi et al.l ()2005r . 
This provides a second method to estimate masses based 
on reverberation mapping. 

We measured the H/3 line flux for the seven NLSls from 
our sample with available optical spectra and used the 
results to compute the black-hole mass in each object 
according to (PJ. The values obtained for L(H/3) and 
the resulting estimates of Mbh are included in Table El 
The average unweighted black-hole mass was found to be 
(log(MB//)) = 6.4, again consistent with the notion that 
black holes in NLSls are statistically less massive than 
those found in broad-line Seyferts. 

4.2. The Mass-Luminosity Relation 

As discussed in § 1, it is generally accepted that the 
slope of the correlation between Mbh and host bulge lu- 
minosity is consistent for both non-active galaxies and 
"typical" Seyfert galaxies. However, it remains a matter 
of debate as to whether the relation can be extended in 
the same form to the NLSl subset. Using the results 
obtained from our profile decompositions, we now derive 
estimates of Mbh under the assumption that the slope 
of the correlation is the same for all Seyfert galaxies, in- 
cluding NLSls. We then compare these mass estimates 
to those obtained in § 4.1. The results are presented 
in this manner in order to avoid a priori assumptions 
about the accuracy of the black-hole masses obtained 
from the mass-radius relation. While reverberation map- 
ping masses are generally secure, only a handful objects 
that could be classified as NLSls have been directly ob- 
served using this technique. It is conceivable that other 
factors such as the orientation of the broad-line region 
could affect the observed line widths, potentially leading 



to underestimates of the black-hole masses. Such issues 
do not affect the mass-bulge luminosity relation. 

The correlation between Mbh and the mass of the 
spheroid component of the host galaxy is frequently 
quantified in terms of the i?-band magnitude (e.g. 
iFerrarese^ Ford 20051) . However, the blue luminosity 
is biased by the presence of young, hot stars and does 
not necessarily best represent the underlying mass dis- 
tribution. Near-IR wavelengths are less sensitive to this 
bias, and are less affected by dust extinction. This is evi- 
dent from a comparison of the mean extinctio n curves in 
B and K' . Using values presented in lMathi3 11990), the 
ratio of B- to if'-band optical depths is approximately 
12.3. Based on the extinction values for our target fields, 
we estimate that the attenuation of light in B due to the 
presence of Galactic dust is more than 11 times greater 
than that in K' . Hence the scatter in the correlation can 
be reduce d by e mploying near-IR data, as first demon- 
strated bv lMarc oni & Hunt (2003). 

We calibrated the mass-luminosity relation for near- 
IR passbands using a sample of 16 bright elliptical 
galaxies with secure mass estimates obtained by resolv- 
ing the black-ho le sphere of influence as reported by 
iFerrarese fc ForJ (^005). The assumption that ellipticals 
are kinematically similar to the bulge components of spi- 
ral galaxies and therefore follow the same correlation is 
2 supported by observational evi dence for consistent Fun- 
^tfoiental Plane relationships (iK ormendy fc Illi ngworthI 
IT981 iFalcon-Barroso. Peletier. fc Balceiy I2002D 
eluding spiral galaxies from the sample eliminates the 
need for separating the total integrated flux into bulge 
and disk components, thus minimizing the scatter in the 
correlation. Furthermore, we reject several ellipticals 
with secure Mbh estimates for the following reasons. 
Cygnus A is known to be a member of a rich cluster of 
galaxies, and thus may be influenced by a complex grav- 
itational potential. Two other galaxies, NGC 821 and 
NGC 2778, have estimated spheres of influence that are 
smaller than the spatial resolution of the data. For each 
remaining object, the total magnitude in the 2MASS 
near-IR bandpasses of J and Ks were extracted from the 
XSC. It is important to note that the response curves 
of the 2MASS filter set are not identical to those of the 
CFHT filters used for our NLSl observations. The nec- 
essary photometric transformations are now discussed. 

Absolute magnitudes were calculated for all but two of 
the ga laxies using distance estimates from Tq^^^ t alJ 
()2001t) . As data were not available in Tonr v et al.l for 
NGC 6251 and NGC 7052, we determined luminosity dis- 
tances using redshifts obtained from the Third Reference 
Catalog of Bright Galaxies (RC3)( dc Vaucouleurs et aO 
Il991[l . A summary of the properties of the galaxies used 
to calibrate the mass-luminosity relation is presented in 
Table 11 

In Figures El and El we plot the logarithm of Mbh, 
in solar masses, versus the absolute magnitude in J and 
Kg for the sample of elliptical galaxies. The best-fit 
slopes to the data were determined from a linear least- 
squ ares fit, weighting: the parameters by their uncertain- 
ties (|Press et al.ll992|) . The mass- luminosity correlations 
in J and Ks have the following forms: 

log Mbh = (-0.43 ±0.04)(Mj-^ 22.52) (8.17 ±0.09) 

(4) 



The Central Engines of NLSls 



7 



\ogMBH = (-0.44±0.04)(Mk +23.45) + (8. 17± 0.09). 

... 

Our results are consistent within unc ertainties with both 
the if^-band fit presented in|Dong fc De Robertis ( 2006), 
o btained using the s ame s ample criteria, and the results 
of iMarconi fc Hiinil ||200.1D . whose sample includes spiral 
galaxies. 

In order to use these correlations, it was necessary to 
convert the bulge magnitudes obtained by GALFIT from 
the Mauna Kea Observatories (MKO) photometric sys- 
tem to the 2MASS system. Color transformations for 
various filter sets are available through the 2MASS web- 
site^. The transformations for the MKO system are: 



Of the 11 NLSls included in our sample, published mor- 
phological classes were available for seven galaxies. Pre- 
liminary classifications were made for two additional 
galaxies, MCG 08.15.009 (T=-l) and PG 1016+336 
(T=l), based on a visual inspection of our high-spatial 
resolution data. Morphological information was not 
available for Mrk 734 and IRAS 04596-2257. Applying 
@ to the sample, bulge magnitudes were calculated and 
compared to those derived by GALFIT. The results are 
included in Table [3 While the scatter in the measure- 
ments for individual galaxies is substantial, the mean 
bulge magnitudes are consistent to within 0.02 magni- 
tudes, further establishing that the values generated by 



K 



s, 2MASS 



GALFIT do not suffer from any systematic errors 
i^MKO + (0.002±0.004) + (0.026±0.006)(J-i4:) MKO 

(6) 5. DISCUSSION 

{J-Ks)2MASs = (1. 037+0. 009)(J-if)M_ft:o+(-0. 001+0. 006). The black-hole masses presented in the previous sec 



(7) 

However, th ese transformations are vali d for the K filter 
and not K' . 'Wains coat &: Cowiel (TM^ compare the K' 
filter to K, and suggest a color transformation of the form 



K' ~K ^ (0.22 + Om){H - K). 



(8) 



tion are plotted as a histogram in Figure [TBI There is a 
clear offset between values obtained using different cor- 
relations. We find that the estimates of Mbh calculated 
from the near-IR bulge luminosities are, on average, more 
than one order of magnitude larger than those obtained 
from the mass-radius relation. The measurement un- 
certainty for each object is typically 0.25 dex, which is 
insufficient to account for the difference in the results. 
Therefore, one (or more) of the assumptions made in 
formulating the correlations must not be applicable to 
this particular subset of active galaxies. In the follow- 
ing subsections, we consider the evidence supporting the 
extension of each correlation to NLSls, and discuss the 
impact of these results on our understanding of these 
objects. 

5.1. Are Black Holes Less Massive In NLSls? 

The most common interpretation of the observed spec- 
tral characteristics of NLSls is that the central black 
holes are less massive than those found in broad-line 
Seyfert galaxies. This hypothesis is supported in the 
literature by studies that employ the mass-radius cor- 
relation as calibrated by reverberation mapping studies 
(e.g. [Grupe & Mathur 2004). Indeed, we find that for 
our sample the mean value of Mbh determined from Q 
and (0) is (log(MBij)) = 6.4. If we accept the virial mass 
estimates obtained from the mass-radius relation, the re- 
sults from our profile decompositions provide evidence 
that the galaxy bulges in our sample are more luminous 
than predicted by the mass-luminosity relation. This is 
illustrated in Figures El and El in which we compare the 
NLSl black- hole masses as determined from the mass- 
luminosity relation and the mass-radius relation. The 
estimates of Mbh derived from the mass-radius relation 
lie well below the best-fit line calculated from ellipticals 
with secure Mbh estimates. 

It has been proposed that the departure of NLSls from 
the mass-luminosity relation is the result of the contri- 
bution of excess flux from eit her bursts of nuclea r star 
formation or AGN emission (INelson et al.ll200l . We 
exercised great care in subtracting the nuclear compo- 
nent from our high-resolution imaging data, however. 
Any residual emission inadvertently included in measur- 
ing the bulge magnitude would be negligible, and cer- 
tainly would not alter the subsequent estimates of Mbh 
by one full order of magnitude. Rather, any errors in the 
^ http://www.ipac.caltcch.edu/2mass/relcascs/allsky/doc/scc6_4b.htiflluclear Subtraction routine would likely involve an over- 



In the absence of //-band data for our sample galax- 
ies, we adopted a mean (H-K) color estimated from 



IFioc fc Rocca-Volmerangd (|l999ft . Using a sample of ap- 
proximately 1000 galaxies, the authors computed near- 
IR colors as a function of host morphology. They deter- 
mine a median (H-K) for ellipticals of approximately 0.20 
+ 0.01, rising to 0.27 + 0.01 for spirals. Assuming that 
the bulge colors for our galaxies are similar to that of 
ellipticals, a conversion from K' to K is possible. Com- 
bining this result with the color transformation between 
K and if,, it can be seen that magnitudes measured in K' 
and Kg are the same within photometric errors. This is 
supported by an inspection of the transmission fu nctions 
for K ' and Kg filter, which are very similar (Ki m et al.l 
I2OOI . A non-trivial transformation between Jmko and 
J2MASS is required of course. 

It was also necessary to adjust our bulge magnitudes by 
applying modest if-corrections to our data. Corrections 
were taken from jPoEgianti (1993) who used evolution- 
ary synthesis models based on redshift and as a function 
of host morphology. When necessary, corrections were 
estimated from a linear interpolation between redshift 
bins of width 0.02. The measured bulge magnitudes were 
then inserted into (^J and ^ to estimate Mbh for each 
NLSl target. The results are presented in Table [31 The 
unweighted mean mass was found to be (log(MB//)) = 
7.9. 

As a consistency check, we applied the 
iSimien k. de VaucouleursI ljl986|) relationship between 
the luminosity bulge-to-disk ratio and Hubble stage. 
This method is less reliable than our more accurate 
component decomposition due to the large dispersion 
in the relation. Using a sample of 118 spiral galaxies, 
representing morphologica l stages SO/ a through Sc, 
iDong fc De Robertis! ()200 6D have calibrated the relation 
for the 2MASS ifg-band. The relation is represented by 
the following cubic equation: 

Am^- = 0.297(r+5)-0.04(r+5)2+0.0035(r+5)^ (9) 



8 



Ryan et al. 



correction, leading to underestimates of the black-hole 
masses. The estimated bulge scale radii are typically 
larger than 0.5", and all are greater than the nominal 
PSF FWHM of 0.2". Therefore, we are confident that 
the derived parameters, including the bulge magnitudes, 
are both physically reasonable and accurate. 

If black holes are less massive in NLSls, there must be 
a physical explanation for the breakdown of the mass- 
luminosity relation. One interpretation is that there is 
an inherent time delay between an increase in the mass 
of the bulge through star formation and the growth of 
the central black hole due to accretion. In this sce- 
nario, sufficient time passes after a gravitational interac- 
tion between galaxies for the bulge to virialize and stars 
to form before gas inflow triggers significant nuclear ac- 
tivity. One could imagine that accretion onto the cen- 
tral black hole is then halted by feedback mechanisms as 
the mass becomes consistent with that of a "normal" 
Se yfert galaxy. Numerical simulations of interactions 
bv lDi Matteo. Springel. &: Hernauis"3 lj2005|) suggest that 
such galactic nuclei undergo a period of rapid accretion 
towards the end of a merger event due to the elevated gas 
content near the core. It is during this phase that the 
central black holes can increase in mass by a significant 
fraction, and the target emerges as an AGN. During this 
period of rapid growth, the accretion rate can approach 
the Eddington limit. If NLSls are indeed evolutionarily 
young objects at the start of the AGN duty cycle, these 
simulations constrain the NLSl phase to a time relatively 
soon after a gravitational interaction, a time when one 
might expect significant morphological asymmetries in 
the host galaxies and particularly the centers. 

In order to search for such morphological evidence, 
we divided the reduced K'-band mosaics by the associ- 
ated J-band images. We then applied an unsharp mask 
by convolving the quotient image with a Gaussian func- 
tion with a standard deviation of 1.5". No evidence was 
found of light asymmetries or tidal features that would 
be indicative of a late-stage merger. Furthermore, we 
note that spectroscopic detections of Seyfert nuclei in 
systems currently undergo ing a stro ng interaction (e.g. 
Liu fc KcnnicuttI [T995t IVanzi. Alonso-Herrero. fc Riekd 
1998.') are difficult to explain if bulge growth and AGN 
fueling are temporally well separated. 

To further investigate the claim that NLSls are evo- 
lutionarily young, we compared the near-IR bulge colors 
of the objects in our study to a subset of non-active, 
early-type galaxies ta ken from the sample examined by 
iDong &: De RobertisI ((2006). Three non-active galax- 
ies were selected for each NLSl, matched according to 
morphological class and absolute magnitude of the total 
galaxy. Of the 27 objects selected, we note that nine are 
actually identified as containing a LINER nucleus. The 
mean bulge color for this subset of galaxies is equivalent 
to the mean for the other galaxies in our non-active sam- 
ple. Therefore, their inclusion does not adversely affect 
our results. Using 2MASS XSC imaging data, the {J- Kg) 
colors of the bulge components of the non-active objects 
were computed with GALFIT. The galaxy parameters 
derived from the fits, including bulge magnitudes, are 
presented in Table El The bulge {J- Kg) colors for the 
non-active galaxies and NLSls are plotted against the 
absolute bulge magnitudes in Kg in Figure lTBI The mean 
bulge color for the non-active sample is {J-Ks) = +1.10 



± 0.04, while for the NLSl sample, the mean color was 
found to be (J-Ks) = +1.85 ± 0.58. 

In the absence of two-dimensional profile decomposi- 
tions for a matching sample of broad-line Seyferts, we 
comp are our NLSl colors to results culled from the liter- 
ature. iHunt et alJ |T997.1 measured the near-IR colors of 
galaxies hosting Seyfert and starburst nuclei. The mean 
color of the inner 3 kpc of the Seyfert 1 galaxies in their 
sample was found to be (J-K) ~ +1. However, the inner 
4" were masked to exclude contamination from the active 
nucleus; hence the measured colors more accurately rep- 
resent the inner disk component rather than the bulge. 
Their use of an early-generation InSb 62x58 IR array 
also limits the comparison to our high-spatial resolution, 
AO-corrected data. 

A comparison of the near-IR colors of Seyfert and 
non-a ctive spiral galaxies is presented in lMarauez et all 
()2000l) . Using 1-dimensional profile decompositions, the 
authors conclude that the bulge and disk components 
of the two galaxy samples are similar in terms of their 
mean {J-K) color and Fundamental Plane relation. How- 
ever there is significant scatter in the measured bulge 
colors of the 18 AGN in their sample, and while the 
mean bulge color was determined to be {J-K) = +1.25, 
several galaxies have a bulge color in excess of +1.5. 
One object in particular, NGC 4253, was found to have 
an exceptionally red bulge, with {J-K) = +4.4. Inter- 
es tingly, while th i s gala xy is classified as a Seyfert 2 
in iMara uez et al J (120001). it i s iden tified as a NLSl in 
'Gren shaw. Kraemer fc Gabell (|2003?). 

We find that of the 11 NLSls in our sample, four have 
a bulge color of {J- Kg) ~ +1 to +1.5, consistent with the 
sample of non-active galaxies and other studies of broad- 
line Seyfert galaxies. However, six of the targets have 
a bulge color of {J-Kg) ~ +1.9, and for one, Mrk 1239, 
we measure {J-Kg) — +3.0. We find no evidence for sys- 
tematic effects that would lead to errors in our measured 
colors; on the contrary, the targets with the reddest bulge 
colors are also the model fits with the lowest residuals, 
as parameterized by the reduced values. 

The near-IR spectrum of Mrk 1239 may provide a clue 
to the source of the unusual bulge colors. A signifi- 
cant bump in nuclear continuum emission, centered on 
2.2 /xm , was detected by iR.odrig^iez-Ardila fc Mazzalavl 
()20n6|) . Given the resolution limit of approximately 
1" for their data, the authors constrain the source of 
this emission to the central 380 pc, which is well out- 
side the seeing disk of our imaging data. A similar 
flux excess has been detected in the NLS ls Mrk 766 
(iRodriguez- Ardila. Gontini. fc Viega^l2005|) and I Zw 1 
(Rudv et al. 2000). Gomparing the spectrum of Mrk 
1239 to that of the NLSl Ark 564, they find the lat- 
ter does not exhibit any excess continuum emission in 
the if-band. Hence the feature appears to be present in 
some, but curiously not all, NLSl galaxies. This may 
explain the large scatter in the {J-Kg) bulge colors of our 
sample. 

Because the observed Jfj-band excess is constrained to 
the central regions of the targets, we closely examined 
the inner bulge colors of our targets. We find that some 
of the NLSls in our sample exhibit a gradient, such that 
the bulge colors become redder with decreasing radius. 
This is consistent with the view that an additional source 
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of fLux is detected at 2.2 /im in the inner regions of some 
of our NLSl targets. Consequently, our Ks-ha.nd bulge 
profiles may be influenced by this component, producing 
the small but systematic residual pattern evident in Fig- 
ure El Fitting the peak with a blackbody temperature 
of 1200 K. lR,odriguez-Ardila fc Mazzalavl 1200^ suggest 
that the source of the excess flux is hot graphite dust 
emission near the active nucleus. This explanation, how- 
ever, cannot account for the spatial extent of the emission 
as demonstrated by the observed color gradients. 

We also note that the estimates of AIbh derived from 
the profiles in Kg are generally larger than those obtained 
from the J-band fits. Given the possible presence of an 
additional source of flux near the centers of some NLSls, 
we consider the mass estimates in J to be more robust, 
and caution against the use of Ks-hand data to measure 
NLSl host bulge luminosities. 

In summary, while we cannot formally exclude the low- 
mass interpretation, we find no evidence from our data 
that NLSls as a group have morphological asymmetries 
or extra stellar components distinct from normal Seyfert 
galaxies that would suggest that these objects contain 
a recently formed AGN. It may well be that there is a 
modest flux excess at 2.2 ^m that affects the near-IR 
bulge colors of some, but not all, NLSl galaxies. 

5.2. Could The Mass- Luminosity Relation Hold For 
NLSls? 

The alternative interpretation of our results is that 
the mass-luminosity relation does indeed hold for NLSls 
and that the masses obtained from the correlation with 
broad-line radius are unreasonably low. It is evident from 
Figures [TSl and ITU that the measured NLSl bulge mag- 
nitudes are similar to those of the non-active galaxies 
used to derive the correlations. The estimates of Mbh 
for the NLSls were obtained by interpolating within a 
region for which the correlations are known to be valid 
for b road-line Seyferts. This result is not in agreement 
with iBotte et a l. ( 2004), who suggest that the i3-band 
luminosities of NLSl bulges are lower than those of 
the broad line counterparts by approximately 1.3 mag- 
nitudes. However, as discussed earlier, near-IR wave- 
lengths should in principle provide a more accurate trace 
of the underlying mass distribution than the 5-band. 

We are confident that our spectroscopic measurements 
do not suffer from any systematic errors. Therefore, 
if we accept the extension of the mass-luminosity rela- 
tion to NLSls, we require an explanation for the narrow 
emission lines that characterize these objects. Although 
strong non-virial motions in the broad-line region are 
known to induce asymmetries in the C iv line profiles in 
AGN spectra, their contribution at optical wavelengths 
is insufficient to significantly affect the H/3 line profile 
ijVeron-Cettv. Veron. fc GoncalvesI l200l|) . We have yet 
to directly resolve the broad line region, and rely on in- 
direct evidence that th e motion of the gas is virialized 
ijFerrarese fc Fo rd 2005). If nuclear activity is triggered 
by galaxy mergers, the gas may be heavily influenced 
by non-gravitational forces, partic ularly at the earliest 
stages of the AGN duty cycle (e.g. 'Barn es fc HernauTsl 
^91). It is possible that the broad- line emitting gas 
in the cores of NLSls does not reflect only the gravita- 
tional potential of the central source. Departures from 
virialized motions would almost certainly broaden emis- 



sion lines, increasing the discrepancy with the mass- 
luminosity rela tion. 

As noted bv iPeterson et al.l l|2000D . there are poten- 
tial explanations for the remarkable spectral properties 
of NLSls other than the low-mass hypothesis. It is pos- 
sible that the orientation of the broad-line region could 
influen ce mass estimates derived using the iKaspi et alJ 
l)2005j) relation. Studies of the central engines of active 
galaxies typically involve a large sample of objects with 
a range of spectral properties, and the assumption that 
the broad-line regions are oriented randomly does not af- 
fect the statistical results. On the other hand, NLSls are 
a small and distinct subset of Seyfert galaxies, and the 
assumption of heterogeneity may not be reasonable. If 
NLSl systems preferentially have a low inclination, val- 
ues of Mbh could be underestimated by a factor of three 
or more. 

In § 14.11 it was mentioned that the model geometry of 
the broad- line region could influence estimates of Mbh- 
An isotropic gas distribution is assumed in many studies 
that employ the mass-radius relation (e.g. ^asoi et all 
12000; Wandcl 2002). The exclusion of a disk component 
in these studies could produce virial black-hole masses 
that are a underestimated by a factor of three. As the 
two effects mentioned here are multiplicative, their net 
effect could lead to values of Mbh that are a full order of 
magnitude too low. Therefore, it is possible that our in- 
congruent mass estimates are the result of the combined 
effects of broad-line geometry and orientation, and the 
mass-luminosity relation could indeed hold for NLSls. 

5.3. Are NLSls Fueled By Merger Events? 

Given the discrepancy between the estimates of Mbh 
revealed in this study, it is apparent that we cannot ap- 
ply both the mass-radius and mass-luminosity relations 
to the NLSl subset in the same form as used for broad- 
line AGN and non-active galaxies. The most common in- 
terpretation is a breakdown in the correlation with bulge 
luminosity. We have exercised a great deal of care in mea- 
suring accurate bulge magnitudes. We also have demon- 
strated that the application of the mass-radius relation to 
NLSls may suffer from systematic effects that could lead 
to underestimated virial masses. Furthermore, we find 
no evidence for features that betray recent tidal interac- 
tions or distinct stellar components that would support 
the interpretation that NLSls are evolutionarily young. 

There may be an alternative view of NLSls that pro- 
vides an explanation for their exceptional properties, 
however. There are some interesting, albeit circumstan- 
tial, lines of evidence that point to a mechanism other 
than galaxy harassment to initiate accretion in these ob- 
jects. It has been suggested that as the universe contin- 
ues to expand and cosmic distances increase the role of 
mergers in defining galaxy morphology will become less 
important compared to secular evolution, the gradual re- 
organization of mass du e to internal processes such a s 
bar-driven instabilities ijKormendv fc Kennicutti 12004^. 
The photometric and kinematic properties of galaxies 
that have undergone secular evolution will be distinct 
and recognizable from those which arise from merger 
events. It is conceivable that NLSls represent such ob- 
jects. 

The first piece of evidence in favor of this hypothesis 
involves the shape of the bulge surface brightness pro- 
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files. Bulges formed via secular evolution are believed 
to be more accurately represented by ex ponential pro- 
files than by the de Vaucouleurs r^/^ law ()Balcells et alJ 
I2OO3,') . When fitting the bulge components for the galax- 
ies in our sample, the Sersic index, n, was allowed to vary 
as a free parameter within the GALFIT minimization 
routine. The average of the optimal indices determined 
by GALFIT was 1.4 in K' and 1.6 in J, and n was never 
larger than 2.4. The NLSl bulges were thus found to be 
more consistent with an exponential profile. 

To be fair, the true nature of bulges remains some- 
what ambiguous. In studying the bulge profiles of SO 
to Sbc galaxies using the NICMOS near-IR camera on 
the Hubble Space Telescope, Balcclls et al. (2003) found 
the average Sersic index to be 1.7 ± 0.7. Ground-based 
near-IR observations of the same sample were best fit 
by indices of 4 to 6. The authors suggest that the dis- 
agreement is caused by nuclear point source blending due 
to the effects of atmospheric seeing in the ground-based 
data. It is possible that this problem was mitigated in 
our ground-based data by the excellent angular resolu- 
tion, combined with the application of our empirical tech- 
nique for removing the nuclear emission. 

Another key piece of evidence supporting secu- 
lar evolution is the possible failure of the Mbh- 
(T.s correlation in NLSls. Num erical simulations by 
iDi Matteo. Soringel. fc HernnuistI ((2()05) demonstrate 
that major galaxy mergers produce strong infiows of gas 
that feed the black hole. The energy emitted by the black 
hole expels gas from the central region, limiting star for- 
mation and regulating further accretion. The connection 
between Mbh and bulge velocity dispersion could there- 
fore be a natural product of galaxy mergers. Conversely, 
the departure from this relationship by NLSls may in- 
dicate that mergers are not largely responsible for initi- 
ating nuclear activity. Geometry may thus play a more 
significant role in the narrow II/3 emission lines. 

If accretion in NLSls is not primarily the result of 
gravitational interactions, there must be an alternative 
mechanism for driving gas towards the nucleus. In the 
secular evolution paradigm, this is accomplished by stel- 
lar bars that produce non-axisymmetric gravitational 
potentials, driving gas inwards. Thus, barred early- 
type galaxies are ideal candida tes for secular evolution 
l|Kormendv fcK ennicutt' "2004*). Interestingly, it was 
reported by [Crenshaw. Kraemer & Gabel ( 2003) that 
large-scale bars appear to be more common in NLSls 
than in broad-line AGN. Of the seven objects in our 
sample for which morphological information was avail- 
able, four are classified as barred galaxies. We did not 
detect any such features in our near-IR data that could 
be included in the decompositions, however. Obser- 
vations of early-typ e barred and unbarred systems by 
lAguerri et al.l pOOSfl demonstrate that the presence of a 
bar does not change the galaxy morphological parame- 
ters such as Sersic index and effective radius. 

Given the arguments presented above, we postulate 
that NLSls may represent AGN in which accretion may 
be regulated by internal processes rather than galaxy 
mergers. While there are outstanding issues that have 
not been addressed in this hypothesis, such as the x-ray 
properties of NLSls and the exact cause of the narrow 
emission features, we feel that further investigation into 



the possibility of secular development is warranted. 

6. CONCLUSIONS 

We have studied a sample of NLSl galaxies in order to 
explore the hypothesis that such objects are evolution- 
arily young AGN possessing relatively low-mass black 
holes. Galaxy profiles were separated into their con- 
stituent components by applying the two-dimensional 
decomposition algorithm GALFIT to high-spatial reso- 
lution, near-IR imaging data. Black-hole masses were 
derived using the correlation between Mbh and host 
bulge luminosity, calibrated in J and K' using a sample 
of bright elliptical galaxies with black-hole masses ob- 
tained by resolving the sphere of influence. The scatter 
in the mass-luminosity relation was thus minimized by 
restricting the calibration sample to objects with secure 
Mbh estimates. Furthermore, the near-IR form of the 
relation is less sensitive to attenuation by dust than at 
optical wavelengths. The unweighted mean black-hole 
mass determined from this relation is, in solar units, 
{\og{AlBH)) = 7.9, in line with typical values for broad- 
line Seyferts. We also derived masses for seven of the 
objects in our sample from the correlation between the 
radius of the broad-line region and the monochromatic 
continuum luminosity at 5100 A under the assumption 
that the emission-line gas is virialized and using optical 
spectra taken from the literature. The mean mass calcu- 
lated from this relation is (log(MBjy)) = 6.4, more con- 
sistent with the premise that the black holes in NLSls are 
less massive than those found in their broad-line coun- 
terparts. 

It is clear that there is a significant discrepancy be- 
tween the results obtained from these correlations, with 
the mass-luminosity relation yielding values of Mbh that 
are on average more than one order of magnitude larger 
than those obtained from the mass-radius relation. For- 
mal uncertainties in the calculated masses are typically 
around 0.25 dex, and thus are insufficient to account 
for the conflicting results. Our findings suggest that (at 
least) one of the correlations cannot be extended to the 
NLSl subset in the same form as used for other galaxies. 
It has been conjectured that NLSls may not follow the 
standard mass-luminosity relation due to contamination 
by circumnuclear star formation. We have argued that 
this extra component would not significantly affect our 
data. Furthermore, we noted that the resolution of our 
imaging data provided excellent spatial over-sampling of 
the surface brightness profiles, producing highly accurate 
bulge luminosities. The resulting Mbh estimates were 
calculated by interpolating the relation within a region 
known to be valid for both non-active galaxies and other 
AGN. 

Numerical simulations of galaxy interactions constrain 
the high-accretion rate period that may define the start 
of the AGN duty cycle to a period shortly after a sig- 
nificant gravitational perturbation of the host galaxy. 
To further investigate this possibility, we compared the 
NLSl bulge (J-Kg) color indices to a matched sample 
of non-active galaxies. We determined the mean bulge 
color for the NLSl sample to be ((J - Ks)) = +1.85 
± 0.58. For the non-active galaxy sample, the mean 
color is ((J - Ks)) = -hi. 10 ± 0.04. Estimates of the 
bulge colors for broad-line Seyfert 1 galaxies were taken 
from the literature. A fair comparison to our sample 
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is complicated by the large scatter in the published re- 
sults and the potential for misidentification of NLSls as 
Seyfert 2 galaxies. The NLSl images were also examined 
for light asymmetries or tidal features that may reflect 
a recent gravitational interaction. We find no evidence 
for any significant morphological asymmetries or extra 
stellar components that would suggest these galaxies are 
evolutionarily young. 

We have, however, identified several lines of cir- 
cumstantial evidence suggesting that nuclear fueling in 
NLSls may be strongly influenced by secular processes 
within the host galaxy. The bulge light profiles obtained 
from our two-dimensional decompositions are all nearly 
exponential in shape, a key characteristic of secular evo- 
lution. NLSls may not follow the standard Mbh-c^s re- 
lation, which has been shown to be a natural product 
of galaxy mergers, while an enhanced fraction of stellar 
bar structures in NLSls may provide a natural source 
for fueling the nucleus without invoking a recent gravita- 
tional interaction. It is thus fair to acknowledge that we 
do not yet completely understand the evolution of this 
important subset of low-redshift active galaxies. 
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TABLE 1 
NLSl Galaxy Parameters 



Target ID 


Distance^ 


Morph. Class. 


Total Magnitude'' 


Optical SMA 


K=+20 Isophotal Radius 




(Mpc) 




J 


K, 


(arcsec) 


(arcsec) 


Mrk 335 


111.1 


so/a 


+12.03 ± 0.04 


+10.06 ± 0.03 


9 


15.1 


Mrk 359 


74.8 


SBOa 


+11.61 ± 0.02 


+10.46 ± 0.03 


18 


16.6 


Mrk 618 


153.6 


SB(s)b pec 


+11.82 ± 0.03 


+10.37 ± 0.03 


27 


18.1 


Mrk 705 


125.7 


SO? 


+11.94 ± 0.03 


+10.74 ± 0.04 


21 


16.1 


Mrk 734 


217.6 


Compact 


+13.04 ± 0.04 


+11.76 ± 0.06 


18 


7.6 


Mrk 1044 


70.7 


SBO 


+11.78 ± 0.03 


+10.47 ± 0.03 


21 


14.0 


Mrk 1126 


45.7 


(R)SB(r)a 


+11.09 ± 0.02 


+10.14 ± 0.04 


54 


26.1 


Mrk 1239 


85.8 


E-SO 


+11.87 ± 0.02 


+9.60 ± 0.02 


12 


10.1 


IRAS 04596-2257 


176.5 




+ 13.21 ± 0.04 


+12.03 ± 0.06 


9.3 


6.4 


MCG 08.15.009 


105.0 




+ 12.42 ± 0.05 


+ 11.28 ± 0.07 


18 


14.1 


PC 1016+336 


105.6 


Compact 


+13.03 ± 0.05 


+11.67 ± 0.06 


9 


11.6 



^Distances were calculated using redshift information obtained from the literature. 

^Total galaxy magnitudes and isophotal radii at if = +20 are taken from the 2MASS XSC; optical radii are taken from the literature. 



TABLE 2 
Summary of Observations 



Target ID Coordinates Exposure Time" Mean Airmass 





R.A. (J2000) 


Decl. (J2000) 


J 


K' 


J 


K' 


Mrk 335 


00 06 19.5 


+20 12 10.5 


1200 


80 


1.18 


1.08 


Mrk 359 


01 27 32.6 


+19 10 43.8 


1200 


480 


1.04 


1.13 


Mrk 618 


04 36 22.2 


-10 22 33.8 


1200 


120 


1.17 


1.16 


Mrk 705 


09 26 03.3 


+12 44 03.6 


1680 


480 


1.04 


1.14 


Mrk 734 


11 21 47.1 


+11 44 18.3 


480 


240 


1.01 


1.02 


Mrk 1044 


02 30 05.4 


-08 59 52.6 


1920 


480 


1.16 


1.27 


Mrk 1126 


23 00 47.8 


-12 55 06.7 


1680 


1680 


1.40 


1.63 


Mrk 1239 


09 52 19.1 


-01 36 43.5 


480 


240 


1.09 


1.09 


IRAS 04596-2257 


05 01 48.6 


-22 53 23.2 


1680 


1920 


1.41 


1.44 


MCG 08.15.009 


07 51 51.9 


+49 48 51.6 


1920 


1920 


1.16 


1.16 


PG 1016+336 


10 19 49.5 


+33 22 03.7 


1680 


480 


1.03 


1.04 



Note. — Units of right ascension are hours, minutes, and seconds; units of declination are 
degrees, arcminutes, and arcseconds. 

^Exposure times represent total on-target integration time. 



TABLE 3 

Photometric and Structural Parameters of Components 



Target ID J-band Results iiTs-band Results 





M(bulge)=' 


M(disk) 


B/D 


Sersic n 


r (kpc)'' 


M(bulge) 


M(disk) 


B/D 


Sersic n 


r (kpc) 


Mrk 335t 


-21.29 


-21.69 


0.70 


1.02 


0.25 


-23.33 


-22.91 


1.43 


1.23 


0.15 


Mrk 359 


-20.42 


-22.38 


0.16 


1.92 


0.23 


-22.10 


-23.97 


0.24 


2.13 


0.23 


Mrk 618 


-21.74 


-23.01 


0.31 


1.20 


0.22 


-23.83 


-27.00 


0.05 


1.55 


0.26 


Mrk 705t 


-22.47 


-24.66 


0.13 


2.40 


0.93 


-23.97 


-25.87 


0.17 


2.28 


0.54 


Mrk 734t 


-22.34 


-25.58 


0.05 


2.10 


0.60 


-23.82 


-24.32 


0.63 


1.47 


0.33 


Mrk 1044 


-20.78 


-21.35 


0.59 


0.96 


0.23 


-22.83 


-22.49 


1.36 


0.83 


0.17 


Mrk 1126 


-20.15 


-22.59 


0.11 


1.64 


0.29 


-21.35 


-22.72 


0.29 


1.75 


0.21 


Mrk 1239 


-21.72 


-22.63 


0.43 


1.65 


0.17 


-24.76 


-25.60 


0.46 


1.08 


0.19 


IRAS 04596-2257 


-21.42 


-21.90 


0.64 


1.16 


0.37 


-23.56 


-23.49 


1.08 


0.81 


0.60 


MCG 08.15.009 


-21.70 






1.45 


1.02 


-22.66 






1.15 


0.56 


PG 1016+336 


-20.41 


-21.06 


0.55 


1.27 


0.19 


-22.32 


-22.05 


1.29 


0.94 


0.12 



Note. — Galaxies with multiple image sets acquired in Ks arc denoted by For these objects, the mean fit parameters are 
listed. 

^Absolute magnitudes of galaxy components derived from GALFIT output. 
'^Bulge scale radius in kpc. 
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TABLE 4 

Galaxies Used to Calibrate Mass-Luminosity Relation. 











A'f 
Ivl ^ 


7 






Mr- 




























IC 1459 


28.8 


26.0 lb 11.0 


-24, 


.62 


± 


0, 


.28 


-25, 


.50 




0.28 


NGC 221 


0.77 


0.025 lb 0.005 


-18, 


.21 


± 


0, 


.08 


-19, 


.36 


± 


0.08 


NGC 3377 


10.9 


i.OO_q;^ 


-21, 


.93 


± 


0, 


.09 


-22, 


.76 


± 


0.09 




LV.O 


i.vjO it U. 1 O 


-zz. 




4_ 
It 


u. 


1 1 
. i i 


-23, 


,81 




nil 
U.li 




22.5 


1-^-0.6 


-22, 


.83 


-|- 


0, 


.14 


-23, 


.67 


± 


0.14 


NGC 4261 


31.1 


5.4 ± 1.2 


-24, 


.28 




0, 


.19 


-25, 


,21 


± 


0.19 


NGC 4291 


25.3 


o 1+0.8 
•^■'^-2.3 


-22, 


.69 


± 


0, 


.32 


-23, 


,62 


± 


0.32 


NGC 4S74 


17.7 


'-' -6.7 


-24, 


.15 


-|- 


0, 


,11 


-25, 


.04 


± 


0.11 


NGC 4473 


15.3 


-I -1+0.5 
1-1-0.8 


-22, 


.91 


± 


0, 


.13 


-23, 


.78 


± 


0.13 


NGC 4486 


15.8 


35.7 ± 10.2 


-24, 


.29 


± 


0, 


,16 


-25, 


,18 


± 


0.16 


NGC 4564 


14.5 


0.56tH^8 


-21, 


.97 


± 


0, 


,17 


-22, 


,89 


± 


0.17 


NGC 4649 


16.4 


20.01^:° 


-24, 


.43 


± 


0, 


.15 


-25, 


,35 


± 


0.15 


NGC 4697 


11.4 


1 7+0.2 
1- '-0.3 


-23, 


.08 


± 


0, 


.14 


-23, 


,94 


± 


0.14 


NGC 5845 


24.7 


2 4+0.4 


-21, 


.94 


± 


0, 


,21 


-22, 


.88 


± 


0.21 


NGC 625lt 


99.4 


5.9 ± 2.0 


-25, 


.01 


± 


0, 


,22 


-25, 


.99 


± 


0.22 


NGC 7052t 


62.5 


n 7+2.6 
' -1.5 


-24, 


.55 


± 


0, 


,22 


-25, 


,45 


± 


0.22 



^Galaxy distances obtained from iTonrv et alj 120011) except for two objects 
denoted by ^ . Distances for these galaxies determined using redshift information 
from the RC3. 
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TABLE 5 
Black Hole Mass Estimates 



Target ID 


FWHM{H/3)=' 


ALa{5100 A) 


L(H/3) 


log A/fli/(Mass- 


-Radius Relation)'' 


log Mbh 


(M, 


ass-Luminosity Relat: 




(km s^-"-) 


(10'''* erg s-i) 


(10''2 erg s-i) 


La(5100 A) 


L(H/3) 


J-band 


Ks 


-band 


T Type 


Mrk 335 


1640 


0.640 ± 0.006 


1.0 ± 0.1 


6.81 


6.99 


7.64 




8.12 


8.37 


Mrk 359 


480 


0.193 ± 0.004 


0.07 ± 0.01 


5.38 


5.28 


7.27 




7.71 


7.82 


Mrk 618 


2300 










7.83 




8.34 


8.31 


Mrk 705 


1990 


0.273 ± 0.004 


0.33 ± 0.04 


6.72 


6.88 


8.15 




8.40 


8.26 


Mrk 734 


1820 


0.728 ± 0.007 


0.6 ± 0.1 


6.94 


6.94 


8.09 




8.33 




Mrk 1044 


1280 


0.153 ± 0.003 


0.19 ± 0.02 


6.17 


6.36 


7.42 




7.90 


7.76 


Mrk 1126 


2500 










7.15 




7.25 


7.49 


Mrk 1239 


910 


0.163 ± 0.004 


0.24 ± 0.04 


5.89 


6.13 


7.83 




8.75 


8.41 


IRAS 04596-2257 


1500 










7.70 




8.22 




MCG 08.15.009 












7.82 




7.82 


7.92 


PG 1016+336 


1600 


0.073 ± 0.001 


0.059 ± 0.006 


6.14 


6.28 


7.26 




7.68 


7.61 


Average 








6.29 


6.41 


7.65 




8.05 


7.99 



^Values for the FWHM of the broad component of H/3 taken from lVeron-Cettv. Veron. &: Goncalvej 120011) . 
^Black-hole masses are presented in log space, in units of solar masses. 
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TABLE 6 

Non-active Galaxy Parameters 



Target ID Distance'' J-band Results ifs-band Results 

(Mpc) M(bulge)'' M(disk) M(bulge) M(disk) 



TP "^f^fi 


1 S 

lo. 


2 
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9Q 
.zy 


9*^ 
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SI 
.O -L 


9^^ 
-zo. 


12 


TP f^9n 




n 

.u 


-22 


.ou 


-zo. 


^9 
.oz 


9*^ 

-ZOi 


.^o 


-24 


22 


NPP 9*^fi 


fi9 
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■ OO 


-24 


.UO 


9*^ 

-zo. 


.Ol 


-24 


.OO 


NPP 91ztfi 


1 7 


2 


-22 


42 


-22 


"1 n 


-24 




-21 


A7 

.4: 1 


IN \_xv_' 1 1 O 


1 7 


n 
u 


-21 


. UO 


9*^ 

-zo. 


1 n 

. -LU 


-22 


9S 
.zo 


9*^ 

-zo. 


Q1 


l\rPP 97S9 
IN V_x v_' z / oz 


O I . 


Q 

o 


-21 


,zo 


-22 




-22 


.ou 


9*^ 
-zo. 


.oz 


NPP 9QS'S 


22, 


4 


-22, 


.81 


-22 


.59 


-23 


.77 


-23 


.41 


NPP "^1 fifi 


22, 


Q 


-21, 


.83 


-22 


.64 


-22 


92 


-23 


.40 


NPP "^1 (SQ 


19. 




-22, 


.74 


-22 


.34 


-23 


.82 


-22 


.93 


NPP "^1 QO 


22, 


4 


-21, 


.76 


-22 


.41 


-22 


.94 


-23 


.47 


NPP '^^f\A 


9fi 


c: 

. o 


-21 




-22 


fin 
.uu 


-22 




9*^ 

-zo. 


44 


NPP '^^R'^ 


O'i. 


n 


9n 

-zu. 




-22 


.OO 


-21 




9*^ 
-zo. 


42 


NPP 'HTO^ 

IN V_r v_j O i UO 


17. 


.0 


-21, 


.36 


-21 


.11 


-22 


.40 


-21 


.88 


NGC 3729 


17. 


,0 


-18, 


.90 


-20, 


.79 


-20, 


.06 


-21, 


.65 


NGC 3884 


91. 


,6 


-22. 


.46 


-24, 


.20 


-23, 


.46 


-25, 


.13 


NGC 3898 


21. 


,9 


-22. 


.47 


-22, 


.28 


-23, 


.43 


-23, 


.10 


NGC 4064 


16. 


,9 


-20. 


.31 


-21, 


.38 


-21, 


.49 


-22, 


.05 


NGC 4192 


16. 


,8 


-20. 


.87 


-23, 


.21 


-22, 


.11 


-24, 


.07 


NGC 4293 


17. 


,0 


-19. 


.59 


-21, 


.73 


-20 


.91 


-22 


.76 


NGC 4369 


21. 


,6 


-20, 


,34 


-21, 


.63 


-21, 


.30 


-22, 


.44 


NGC 4457 


17. 


,4 


-21, 


,79 


-21, 


.97 


-22 


.80 


-22 


.68 


NGC 4643 


25. 


,7 


-23, 


,58 


-22 


.11 


-24. 


.49 


-23 


.07 


NGC 4665 


17, 


,9 


-21, 


,13 


-22, 


.67 


-22 


.07 


-23 


.44 


NGC 5377 


31. 


,0 


-21, 


,95 


-22, 


.85 


-22, 


.96 


-23 


.68 


NGC 5701 


26. 


,1 


-23, 


.25 


-20, 


.66 


-24, 


.08 


-21 


.97 


NGC 6340 


22. 


,0 


-20, 


.78 


-22, 


.21 


-21, 


.72 


-23, 


.03 


NGC 6654 


29. 


,5 


-21, 


.41 


-22, 


.60 


-22, 


.37 


-23, 


.47 



^Distances were calculated using redshift information obtained from the 
literature. 

'^Absolute magnitudes of components derived from GALFIT model fits to 
2MASS data. 
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Fig. 1. — J-band surface brightness profile for Mrk 1239. The top panel shows the relative contribution of each component. In the lower 
panel, the residuals (observation minus model) are plotted. 
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Fig. 2. — J-band surface brightness profile for Mrk 335. 
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Fig. 3. — J-band surface brightness profile for Mrk 1126. 
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Fig. 4. — J-band surface brightness profile for Mrk 618. 
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Fig. 5. — J-band surface brightness profile for Mrk 359. 
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Fig. 6. — J-band surface brightness profile for Mrk 1044. 
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Fig. 7. — J-band surface brightness profile for Mrk 705. 
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Fig. 8.— J-band surface brightness profile for MCG 08.15.009. 
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Fig. 9.— J-band surface brightness profile for PG 1016+336. 
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Fig. 10. — 7-band surface brightness profile for Mrk 734. 
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Fig. 11.— 7-band surface brightness profile for IRAS 04596-2257. 
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Fig. 12. — _ffs-band surface brightness profile for Mrk 1126, illustrating the residual pattern seen near the centers of many of the model 
fits in Kg. 
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Fig. 13. — Correlation between black hole mass and host galaxy bulge luminosity, calibrated for ./-band data. Open circles (colored blue 
in electronic edition) denote elliptical galaxies with secure estimates of Af/jj/. NLSl mass estimates computed using this correlation are 
plotted as (red) squares. NLSl black-hole masses derived using the mass-radius relation are denoted by (green) triangles. The (blue) solid 
line represents the weighted linear least-squares fit to the elliptical galaxy data. 
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Fig. 14. — Correlation between black hole mass and host bulge luminosity, calibrated for if^-band data. 
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Fig. 15. — Histogram showing the distribution of black hole masses derived in this study. 
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Fig. 16. — Comparison of the bulge {J-Kg) colors for a sample of non-active galaxies and NLSls. 
galaxy sample. NLSl galaxies are denoted by filled triangles. 



Open circles denote the non-active 



